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ABSTRACT 

Aims. To calculate the power budget for electron acceleration and the efficiency of the plasma emission mechanism in a post-flare 
decimetric continuum source. 

Methods. We have imaged a high brightness temperature (~ 10 9 K) post-flare source at 1060 MHz with the Giant Metrewave Radio 
Telescope (GMRT). We use information from these images and the dynamic spectrum from the Hiraiso spectrograph together with 
the theoretical method described in Subramanian & Becker (2006) to calculate the power input to the electron acceleration process. 
The method assumes that the electrons are accelerated via a second-order Fermi acceleration mechanism. 

Results. We find that the power input to the nonthermal electrons is in the range 3 x 10 25 -10 26 erg/s. The efficiency of the overall 
plasma emission process starting from electron acceleration and culminating in the observed emission could range from 2.87 x 10~ 9 
to 2.38 x 10 \ 

Key words. Acceleration of particles, Plasmas, Sun: radio radiation, Sun: corona 



1. Introduction 

We observed the 1060 MHz counterpart of an M2.8 flare with 
the Giant Metrewave Radio Telescope (GMRT) on Nov 17 2001. 
This long-duration flare started at 04:49 UT, peaked at 05:25 
UT and ended at 06: 1 1 UT. This flare was accompanied by a 
fast halo coronal mass ejection (CME). The observations of this 
flare have been described in detail in Subramanian et al (2003). 
They studied the initiation of the flare-CME event and found ev- 
idence for the formation of a current sheet below the erupting 
arcade that constitutes the CME. We concentrate our attention 
on a relatively simple morphology, high brightness temperature 
source that dominates the 1060 MHz emission after 05:30 UT. 
This source is part of a decimetric continuum feature that ex- 
tends from ~ 700-1100 MHz. The high brightness temperature 
of this source leads us to postulate a plasma emission mechanism 
for this source. We estimate the power input to the accelerated 
electrons that cause this emission and evaluate the overall effi- 
ciency of the process. 

The plasma emission process is thought to be operative in 
several types of meter wavelength emission such as type 1/noise 
storm emission, type II and type III emission (e.g., Melrose 
1975; Gary & Hurford 2004). The overall emission process 
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comprises of several steps: firstly, it is usually assumed that 
a distribution of nonthermal electrons exists. If these nonther- 
mal electrons have a distribution that is anisotropic in velocity 
and/or physical space, they are capable of generating a high in- 
tensity population of Langmuir waves. These Langmuir waves, 
in turn, coalesce with a suitable population of low frequency 
waves to produce the observed radio emission. The efficiency of 
the Langmuir wave-observed radio emission part of this over- 
all process has been well studied (e.g., Robinson, Cairns & 
Willes 1994; Mitchell, Cairns & Robinson 2003). In the context 
of the plasma emission process, the electron acceleration pro- 
cess that produces the nonthermal electron population in the first 
place has received attention relatively recently (Subramanian & 
Becker 2004; 2006). Subramanian & Becker (2004, 2006) have 
calculated the power input to the electron acceleration process 
and the efficiency of the overall process (starting from electron 
acceleration and culminating in the observed radiation) for noise 
storm continua. We use much more reliable observational inputs 
to calculate such an efficiency for the post-flare decimetric con- 
tinuum we have observed. Accurate estimates of the efficiency of 
the overall plasma emission process are important, for instance, 
in deducing the shock strength from the intensity of observed 
type II emission. If the observed intensity can be well related 
to the energy in the nonthermal electrons, one can place reli- 
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able constraints on the parameters of the shock that produced 
the nonthermal electrons. Apart from the conventional type II 
emission, it is conceivable that meter wavelength emission from 
standing shocks in the vicinity of a reconnection region (Aurass 
et al. 2002) can also arise from the plasma emission process. 

2. Summary of radio observations 

Subramanian et al (2003) show that the 1060 MHz emission 
from this flare arises from two distinct sources in the GMRT 
images. The first two peaks of the 1060 MHz total lightcurve 
(top panel of figure 1) arise from the N-S oriented source shown 
in the first two panels of figure 9, Subramanian et al. (2003). We 
concentrate here on the post-flare 1060 MHz emission (after ~ 
05:30 UT) that corresponds to the third broad peak of the total 
lightcurve (top panel of figure 1). The emission after 05:30 UT 
arises primarily from the southern (S) source shown in the last 
panel of figure 9 of Subramanian et al (2003). This is evident 
from the bottom panel of figure 1, where we have plotted the 
emission arising only from this S source. The 1060 MHz con- 
tours of the S source are shown superposed on a soft X-ray SXI 
image in figure 2. 

Figure 3 shows the dynamic spectrum corresponding to this 
event from the Hiraiso spectrograph. The 1060 MHz emission 
from the southern source imaged with the GMRT (figure 2) is 
part of the decimetric continuum shown in figure 3. The deci- 
metric continuum emission spans a frequency range of approxi- 
mately 700 to 1 100 MHz. 

3. 1060 MHz source size and brightness 
temperature 

We have fitted a two-dimensional Gaussian to the southern 
source shown in figure 2 and derived a full width at half- 
maximum (FWHM) size of 36 x 95 . The beam size in the 
images is 30", and we use this to get a deconvolved source size 
of 20 x 90 . We use the 1 GHz Nobeyama Radio Polarimeter 
(NoRP) lightcurve to determine that the peak post-flare flux 
(which, as we know, arises almost exclusively from the S source) 
is around 430 SFU. The brightness temperature T/, (in K) of a 
source is given by 



T h = 1.22 X 10 



10 



f 2 



(1) 



where S is the flux in sfu, / is the observing frequency in GHz, 
and 9 and (p are the source dimensions in arcseconds. We use 
/ = 1.06, S = 430, 9 = 20 and cf> = 90 to get a peak brightness 
temperature of Tb = 2.7 x 10 9 K for the S source. The value 
for the peak flux S is taken from the 1 GHz Nobeyama Radio 
Polarimeter (NoRP) data, since the GMRT data does not have 
absolute amplitude calibration. The high value of the brightness 
temperature strongly suggests a plasma emission hypothesis for 
the observed radiation. 



4. Power input to electron acceleration process 

As mentioned earlier, it is commonly accepted that the starting 
point for the plasma emission process is a distribution of non- 
thermal electrons that is anisotropic in velocity space, and/or in 
physical space. An example of of the former is the "gap" distri- 
bution, which comprises of a thermal distribution together with a 
nonthermal hump that is centered at a velocity v c that is consider- 
ably greater than v e , the thermal velocity. The name derives from 



the relatively empty part of the distribution between the thermal 
and the nonthermal populations. Subramanian & Becker (2006) 
have invoked a second-order Fermi acceleration process to ac- 
count for the nonthermal electrons. They have self-consistently 
combined the requirements of the acceleration process with the 
requirements on the number of nonthermal electrons (relative 
to the number of thermal ones) to arrive at a upper limit for 
v c . The acceleration timescale is lesser than the Coulomb loss 
timescale for electrons with velocities greater than v c ; these elec- 
trons therefore experience net acceleration, while those with ve- 
locities below v c experience (net) collisional losses and eventu- 
ally join the thermal population. Subramanian & Becker (2006) 
express their results in terms of the dimensionless energy f c that 
is related to v c in the following manner: 



{ c = (l/2)m e v 2 c /kT , 



(2) 



where the quantity m e is the electron mass, k the Boltzmann con- 
stant and T the coronal temperature. The upper limit on £ c is 
denoted by £ max , and is self-consistently determined as a func- 
tion of the spectral index 0:2 of the nonthermal electrons from 
the following equation (Subramanian & Becker 2006), where the 
plasma frequency co p can be taken to be equal to the observing 
frequency for our purposes: 



4.46 



2?™ x 2/2 (3 + a 2 )r(^ ; f milx ) 



'(f max + 1.5) 



V? (3+2a 2 ) 
vl/2 



(3) 



where Erfc is the complementary error function. Under the 
plasma emission hypothesis, continuum emission at the funda- 
mental will be polarized in the sense of the 'o' mode, while emis- 
sion at the second harmonic will be polarized in the sense of the 
'x' mode. Since observed type 1 continua are typically polarized 
in the sense of the 'o' mode, the continuum is assumed to be fun- 
damental emission (Melrose 1980). We follow the same logic in 
assuming that the decimetric continuum radiation discussed in 
this paper is fundamental emission via the plasma emission pro- 
cess. 

4.1. Power input to nonthermal electrons L in 

Once £ max is determined, the minimum nonthermal electron en- 
ergy density U* (ergcrrT 3 ) can be determined using the follow- 
ing equation: 



u. 
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The minimum value of the power input L; n (erg s ') to the 
electron acceleration process is related to U» by (Subramanian 
& Becker 2006) 



L in = 8VD t/», 
where 



1.2An e 

£3/2 T 3/2 ' 
b max 1 e 



(5) 



(6) 



The quantity V is the volume of the acceleration region, and Do 
is a constant with units of inverse time, that is equal to 1/4 of the 
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acceleration timescale (Subramanian, Becker & Kazanas 1996). 
The quantity A is the usual Coulomb logarithm, n e the ambient 
electron density and T e the electron temperature. 

We assume the acceleration volume to be the same as that 
from which the observed emission emanates. We envisage this 
region to be a rectangular column. The dimensions of the face 
of this column would be equal to the observed source size, 
20" x 90", which is equivalent to ~ 15467 km x 69600 km. 
Its depth can be calculated from the bandwidth of the emission 
as evident from the Hiraiso dynamic spectrum (Fig 3). The dy- 
namic spectrum shows that the post-flare emission feature start- 
ing at ~ 05:30 UT that we have imaged with the GMRT at 1060 
MHz ranges from ~ 700 to 1 100 MHz. Using the height-density 
model developed in Aschwanden & Benz (1995), we find that 
this frequency extent corresponds to a depth of ~ 6289 km for 
the emission column. The volume of the acceleration region is 
thus taken to be 

V = 15467 x 69600 x 6289 = 6.77 x 10 12 km 3 . (7) 

We use A = 29, T e = 10 6 K, co p = In x 1060 MHz, an 
ambient electron density corresponding to a plasma frequency 
of 1060 MHz, n e = 1.393 x 10 10 cm 3 . Using these values in 
equations |3]-|7] we get 

10 26 > L m > 3 x 10 25 ergs" 1 (8) 

for -10 < (22 < -5. The power law index a-i is defined in such 
a manner that it has to be < -5 in order for the power in the 
accelerated electron population to remain finite (Subramanian 
& Becker 2004; 2006). Save for values of ai approaching the 
limiting value of -5, L m is not a very sensitive function of az. It 
may be noted that we have assumed the parameters relevant to 
the 1060 MHz layer to hold throughout the emission volume in 
calculating L m . 

5. Power in observed radiation L out 

We now turn our attention to estimating the power in the ob- 
served decimetric continuum. As noted earlier, the 1060 Mhz 
post-flare radiation imaged with the GMRT is part of the deci- 
metric continuum feature starting at around 05:30 UT and span- 
ning a frequency range of ~ 700 to 1100 MHz (figure 3). The 
peak intensity of this post-flare emission is S pea k = 430 SFU 
(figure 1). We assume that this is the intensity of radiation emit- 
ted throughout the bandwidth of 1 100 - 700 = 400 MHz. The ra- 
diation flux is therefore F out = 400xl0 6 x430 SFU = 1. 72x10-" 
WnT 2 . The power in the observed radiation L out can be related 
to F out using (e.g., Elgaroy 1977; Subramanian & Becker 2006) 

iou, = F out R 2 Qe- T , (9) 

where R is the sun-earth distance, Q is the solid angle through 
which the radiation is beamed and t is the optical depth of the 
source. Adopting a value of 0.6 steradian for Q. and 1 for r, we 
get L out = 8.54 x 10 17 ergs" 1 . The efficiency 77 = L oM /L m of the 
overall process is thus 

2.38 x 10~ 8 > 77 > 2.87 x 10~ 9 (10) 

6. Conclusions 

We have imaged a high brightness temperature decimetric con- 
tinuum source at 1060 MHz with the GMRT. The emission from 
this source followed a M2.8 flare on Nov 17 2001. Its bright- 
ness temperature is 2.7 x 10 9 K and its observed dimensions 



are 20 x 90 . The emission is part of a decimetric continuum 
feature spanning a frequency range of ~ 700-1100 MHz. The 
high brightness temperature suggests that this emission is due 
to nonthermal electrons emitting plasma waves that are subse- 
quently converted into the observed electromagnetic radiation. 
We follow the formalism of Subramanian & Becker (2006) in 
calculating the power input L; n to the second-order Fermi ac- 
celeration process that produces the nonthermal electrons. We 
use the source dimensions measured from the GMRT images 
and estimate the column depth of the acceleration region from 
the frequency extent of the decimetric continuum to find that 
10 26 > L m > 3 x 10 25 erg s" 1 . We find that the efficiency (77) of 
the overall process starting from electron acceleration and cul- 
minating in the observed radiation is in the range 2.38 x 10~ 8 > 
77 > 2.87 x 10~ 9 . 

The numbers quoted here for L; n and 77 are the best estimates 
of these quantities for the plasma emission process that we are 
aware of. There have been attempts to calculate these quantities 
for noise storm continua (e.g., Subramanian & Becker 2006), 
but the observational inputs for these estimates were not as good 
as the ones used here. In particular, they used oft-quoted, but 
nonetheless rather rough estimates for the emission volume V 
and the output power L out . These quantities are determined much 
more reliably from observations for the calculations presented 
here. 

Before concluding, it is appropriate to reiterate the follow- 
ing caveats: one of the key observational inputs to this calcula- 
tion is the size of the post-flare decimetric continuum source as 
determined from the 1060 MHz GMRT images. The 20" X 90" 
source is significantly more elongated in one dimension as com- 
pared to the other, suggesting that it might be oriented along a 
loop. The actual transverse dimension might then be somewhat 
smaller than the measured size of 20 . Scattering due to coronal 
turbulence is known to be anisotropic (e.g., Chandran & Backer 
2002). The contrast between the longitudinal and transverse di- 
mensions might thus also be accentuated by preferential scatter- 
ing along the longitudinal dimension. The other important ob- 
servational input we use in our calculations is the column depth 
of the radiating source, which we calculate in an approximate 
manner from the extent of the decimetric continuum feature on 
the Hiraiso dynamic spectrograph (figure 3). In calculating the 
power in the observed radio emission L out , we also assume that 
the flux at 1 GHz is the same as that emitted throughout the ~ 
700-1 100 MHz extent of the decimetric continuum emission. In 
the absence of calibrated multifrequency measurements, these 
are the most reasonable assumptions we can make. The formal- 
ism of Subramanian & Becker (2006) that we have used to calcu- 
late the power fed to the nonthermal electrons, L m accounts only 
for anisotropy in velocity space of the nonthermal electrons (the 
gap distribution). In reality, the nonthermal electrons are likely 
to be anisotropic both in physical space (as in a losscone distri- 
bution) as well as in velocity space (e.g., Thejappa 1991; Zaitsev 
et al. 1997). Finally, we have asumed a coronal temperature of 1 
MK in our calculations. There are indications that the tempera- 
ture might be as high as 2 MK, especially above active regions 
(e.g., Hayashi et al. 2006). As shown in Subramanian & Becker 
(2006), doubling the value of the ambient coronal temperature 
will result in the efficiency estimate of the overall process being 
enhanced by a factor of ~ 10. 

We appreciate insightful comments from the anonymous ref- 
eree that has improved the content of this paper. 
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Fig. 1. 1060 MHz lightcurves from the GMRT images for the M2.8 flare on Nov 17 2001 . Since the GMRT images are not amplitude 
calibrated, we have normalized the amplitudes to the peak of the 1000 MHz lightcurve from the Nobeyama Radio Polarimeter 
(NoRP). Top panel: Total lightcurve for all the sources (figure 8, Subramanian et al. 2003). Bottom panel: Lightcurve for the 
southern (S) source only. It is evident that the high brightness temperature S source is the dominant contributor to the part of the 
lightcurve after ~ 5:30 UT. 
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Fig. 2. Contours for the 1060 MHz southern (S) source overlaid on an SXI soft X-ray image. 




Fig. 3. The dynamic spectrum for this event from the Hiraiso spectrograph. The decimetric continuum emission feature is labelled 
in the figure. 



